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Abstract 

In this paper we study a facility location problem in the plane in which a single 
point (facility) and a rapid transit line (highway) arc simultaneously located in order 
to minimize the total travel time from the clients to the facility, using the L\ or 
Manhattan metric. The rapid transit line is given by a segment with any length and 
orientation, and is an alternative transportation line that can be used by the clients 
to reduce their travel time to the facility. We study the variant of the problem in 
which clients can enter and exit the highway at any point. We provide an 0(n 3 )-time 
algorithm that solves this variant, where n is the number of clients. We also present 
a detailed characterization of the solutions, which depends on the speed given in the 
highway. 



1 Introduction 

Suppose we are given a set of clients represented as a set of points in the plane, and a 
service facility represented as a point to which all clients have to move. Every client can 
reach the facility directly, or use an alternative rapid transit line called highway in order 
to reduce the travel time. The highway is a straight line segment of arbitrary orientation. 
If a client moves directly to the facility, it moves at unit speed and the distance traveled is 
the Manhattan or L\ distance to the facility. In the case where a client uses the highway, 
it travels the L\ distance at unit speed to one point of the highway, traverses with a speed 
v > 1 the Euclidean distance to other highway point, and finally travels the L\ distance 
from that point to the facility at unit speed. All clients traverse the highway at the same 
speed. The highway is used by a client point whenever it saves time to reach the facility. 
Given the set of points representing the clients, the facility location problem consists in 
determining at the same time the facility point and the highway in order to minimize the 
total weighted travel time from the clients to the facility. The weighted travel time of a 
client is its travel time multiplied by a weight representing the intensity of its demand. 

Recent papers have dealt with geometric problems considering travelling distances as a 
combination of planar and network distances. Carrizosa and Rodn'guez-Chfa [TU] intro- 
duced the p-facility min-sum location problem on the plane with a metric induced by a 
gauge and a finite set of rapid transit lines giving the network distance. This problem was 
further developed by Gugat and Pfeiffer [15] , and Pfeiffer and Klamroth [21] . Brimberg et 
al. [TJ.E] studied the location of a new single facility considering given regions of distinct 
distance measures. All these papers consider the well-known Weber problem under a new 



'Departamento Matematica Aplicada II, Universidad de Sevilla, Espafia. {dbanez, iventura}@us . es 

^Universitat Politecnica de Catalunya (UPC), Barcelona, Espafia. mkormanc@ulb.ac.be 

*Escuela de Ingenieria Civil en Informatica, Universidad de Valparaiso, Chile, pablo.perez@uv.cl 



1 



metric. Other papers can be found in the context of the combination of the L\ distance 
with the network distance. Abellanas et al. pQ introduced the time metric model: Given 
an underlying metric, the user can travel at speed v{h) when moving along a highway 
h or unit speed elsewhere. The particular case in which the underlying metric is the L\ 
metric and all highways are axis-parallel segments of the same speed, is called the city 
metric [3]. In the scenario of setting up an optimal distance network that minimizes the 
maximum travel time among a set of points, several problems have been recently investi- 
gated in detail [3,5,9]. Other similar and more general models were studied by Korman 
and Tokuyama [18] . See [T7] and [13] for surveys on highway location and extensive facility 
location problems, respectively. 

A similar problem consisting in simultaneously locating a service facility point and a high- 
speed line (i.e. highway) of fixed length was recently studied by Espejo and Rodn'guez- 
Cln'a [13] and Dfaz-Bahez et al. |12j . The authors considered that highway is a turn- 
pike [18], that is, clients can enter and exit the highway only at the endpoints. A first 
solution was introduced by Espejo and Rodrfguez-Chfa, and after that an improved so- 
lution was given by Dfaz-Bahez et al. The problem aims to minimize the total weighted 
travel time from the demand points to the facility service, and can be solved in 0(n 3 ) 
time [12], where n is the number of clients. Dfaz-Bahez et al. [IT] continued the study 
of this variant by considering the min-max optimization criterion. They minimize the 
maximum time distance from the clients to the facility point. 

In this paper we study a related problem in which the length of the highway is variable, 
that is, it is not fixed in advance as part of the input of the problem, and clients can 
enter and exit the highway at any point. Due to the latter condition, highway is called 
freeway [18]. Since both entering and leaving the highway are allowed at all its points, 
then the structure (i.e. highway) is continuously integrated in the plane. We minimize the 
total weighted transportation time from the demand points to the facility. The problem 
of locating a min-max freeway of fixed length was solved by Dfaz-Bahez et al. |11] in 
0{n log n) time. 

The following notation is introduced in order to formulate the problem. Let S be the set 
of n demand points, / be the service facility point, h be the highway, and v > 1 be the 
speed in which demand points move along h. Given a demand point p, w p > denotes the 
weight of p. The travel time between a demand point p and the service facility /, denoted 
by d h (p,f), is equal to: 



The problem can be formulated as follows: 

The Freeway and Facility Location problem (FFL-problem) Given a set S of n 
demand points, the weight w p > of each point p of S, and fixed speed v > 1, locate the 
facility point / and the highway h in order to minimize the next function: 



Our results. We first show that there exist optimal solutions of the FFL-problem in 
which the highway h has infinite length and the facility point / is located on h. We then 
consider only optimal solutions satisfying these properties. We second show that for all 
demand points p the shortest path from p to the facility point / has one of three possible 
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shapes: (a) p moves directly to /, (b) p first moves vertically to reach h and after that 
moves along h to reach /, and (c) p first moves horizontally to reach h and after that 
moves along h to reach /. For each demand point, the shape of its shortest time path 
to / depends on both the speed v in which demand points moves along h and the slope 
of h. This discretization on the shortest path shapes allows us to simplify the expression 
of df l (p, f) and then to obtain a clear expression of the objective function <&(/, h). Using 
geometric obervations, we reduce the search space of the optimal solutions. This is done 
by considering the grid G defined by all axis-parallel lines passing through the demand 
points. We prove the existence of optimal solutions (/, h) that satisfy one of the following 
two properties: (1) h passes through a demand point and / belongs to a line of grid G, 
and (2) / is a vertex of G. The discretization of the search space permits us to obtain the 
main result of this paper, a general 0(n 3 )-time algorithm to solve the FFL-problem. Our 
algorithm divides the search into two cases that correspond to the above two properties. 
As a surprising result, we prove that when speed v is greater than ~ 1.060660172 the 
algorithm can avoid the search of optimal solutions satisfying property (2) because in that 
case there always exists an optimal solution which holds property (1). This result simplifies 
the algorithm when speed exceeds that bound. We finally present three examples, two of 
them showing that when speed is increased and we keep the same configuration of demand 
points the shapes of the shortest time paths can change. A third example shows that when 
speed is less than ^j^, there exist configurations in which the optimal solution satisfies 
property (2). 

Outline. The discretization on the shapes of the shortest paths from the demand points 
to the facility is stated in Section [2j In Section [3] we show how the search space of optimal 
solutions can be reduced. In Section[3]the algorithm to solve the FFL-problem is presented 
and in Section [5] we give the refinement of it. In Section [6l the examples are presented. 
Finally, in Section we present the conclusions and further research. 

2 Discretization of the shortest paths 

Any solution to our problem will be encoded by a pair of elements (/, h), where / is the 
facility point and h is the highway. Given / and h, we say that a demand point p does 
not use h (or goes directly to /) if dh(p, f) is equal to \\p — f\\±. Otherwise we say that p 
uses h. Given any point u of the plane, let x u and y u denote the x— and y— coordinates 
of u, respectively. 

Claim 2.1 Let p± and pi be demand points using the highway h such that they move in 
contrary directions along h. Let segments s±,S2 Q h denote the portions of h traversed by 
pi and p2, respectively. Segments s\ and S2 have disjoint interiors. 
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Proposition 2.2 There exists an optimal solution of the FFL-problem in which the facil- 
ity point is located on the highway. 



Proof. Let (/, h) denote an optimal solution of the FFL-problem and suppose that / does 
not belong to h. Let h! be a translation of h such that / belongs to U . We select h! so 
that to satisfy a condition that will be stated later. Let p be a demand point. If p does 
not use h then: 

d h >(pj)<\\p-f\\i = d h (p,f). (3) 
Otherwise, if p uses h, let q\,q2 G h be points such that 

dh(p, f) = \\P ~ 91 111 + — g2 ^ 2 + ||(?2 - /Ill- 
Let §3 be the point / + (gi — q%), which belongs to the line containing h' . Observe from 
Claim 12.11 that we can select h' so that point q$ belongs to h! for every demand point p 
using h. Then, by using the triangular inequality with the L\ metric, we obtain: 

dyipj) < \\ P -g 3 ||i + 1,93 ~ 71,2 



= ||p~gl + (g2-/)||l+ " gl ^ 
s- II || 1 1 91 _ 92||2 I, ,,, 

< ||p-gi||H H|g2-/||l 

= d h (p,f) (4) 

From equations Q and (@D we have $(/, h!) = YlpeS w p ' d h>(p, f) < Y, P eS w p ' d h(p, f) = 
<&(/, h). Then the pair (/, h') must be an optimal solution and the result thus follows. □ 



Results similar to Propostion l2.2| stating that the facility point belongs to the correspond- 
ing highway, can be found in |lipi4] . Observe from equations ([2]) and ([7]) that there always 
exists an optimal solution (/, h) to the FFL-problem in which the length of h is infinite. 
We then assume from this point forward that every solution satisfies that the highway is a 
straight line and the facility point belongs to the highway. Observe that this assumption 
does not have negative consequences, due to the fact that in practice a highway of infinite 
length is not possible. In fact, if a solution (/, h) to the problem is such that h is a straight 
line, then (f,h") is also an optimal solution, where h" C h is the segment of minimum 
length such that every demand point both enter and exit h on a point of h" . 

Let a always denote the non-negative angle of the highway with respect to the positive 
direction of the x-axis. Unless otherwise specified, we assume < a < j. Observe that if 
a > I we can, by properties of L\ and L2 metrics, modify the coordinate system so that 
angle a satisfies < a < j. 

Given highway h and a demand point p, let p' be the intersection point between h and 
the vertical line passing through p. Similarly, let p" be the intersection point between h 
and the horizontal line passing through p. Let h p i denote the half-line contained in h that 
emanates from p' and does not contain p", and h p n denote the half- line contained in h 
that emanates from p" and does not contain p' . Notice from the assumption < a < j 
that given h and a demand point p, p' is the nearest point to p on h under the L\ metric. 

The next lemma characterizes the way in which demand points move optimally to the 
facility. Let ip v = ^ — arcsin ( ^ J . Since v > 1 we have < ip v < \ . 
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Lemma 2.3 Given the highway h and the facility point f located on h, the following holds 
for all demand points p. If < a < cp v then p moves first to p' and after that moves to f 
using h. Otherwise, if <p v < a < j, the next statements are true: 

(a) If f £ h p i then p moves first to p' and after that moves to f using h. 

(b) If f £ hp/i then p moves first to p" and after that moves to f using h. 

(c) If f E h \ (hp' U hp") then p moves directly to f without using h. 

Proof. Let p be a demand point and assume w.l.o.g. that p is below h. Consider the 
function giu) := \\p — u\\\ + — J^ 2 for all u € h. Notice that g is convex because it is 
a sum of two convex functions. Then any local minimum of g is a global minimum. Let 
= 1- a. 

Suppose < a < cp v . Given e > small enough, let u\ G h \ h p > and U2 G h p i be the 
points such that \\p — u±\\i = \\p — i*2||l = \\p — p'\\i + V%£- Refer to Fig. [TJ 




Figure 1: Proof of Lemma [2.31 The boundary of the square represented with solid lines is the 
set of points u such that \\p — u\\\ = \\p — p'\\i, and the perimeter of the square represented with 
dotted lines is the set of points u such that \\p — u\\\ = \\p — p'\\i + \/2e. 



Then we have the following: 
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From equations ([5]) and © we conclude that g(p') is the minimum of g. Therefore we 
have dh{p, f) = g(p') and the first part of the lemma thus follows. 

Suppose now ip v < a < \. Then we have three cases: 

Case 1: f G h p >. On one hand we have g[u) > g(p') for all points u G h\h p i . On the other 
hand, if e > is small enough and 112 G h p > is the point such that \\p— U2II1 = \\p— p'Hi+v^e, 
then g(u 2 ) - g(p') = e (V2 - ^j^) > 0. Therefore, d h (p, f) = g(p') and statement (a) 
follows. 

Case 2: f G h p //. Let e > be a small enough value. On one hand, if u± G h p // is the point 
such that ||p = ||p-p"||i + V2e, then g(u{) - g(p") = e (y/2 - > 0. On the 

other hand, if 112 G h \ h p n is the point such that ||p — 1 1 1 = ||j> — — V^E; then 



g{u2) - g{p") 



1 



V2 



smi 



> e 







8111 



<Pv) ■ 



V2 



Therefore, dh(p, f) = g{p") and statement (b) follows. 

Case 3: f G h \ (h p i U h p n). If 9 = then / is one of the nearest points to p on h by 
considering the L\ metric. Thus g(u) = \\p — + ^ U ~J^ 2 > \\p — u\\\ > \\p — /||i = #(/)• 
Otherwise, if 9 > 0, we proceed as follows. On one hand we have g(u) > g(f) for all points 
u G h to the left of /. On the other hand, if e > is small enough and 112 G h is the nearest 
point to / satisfying ||p - u 2 \\i = \\p - f\\ - y/2e, then g(u 2 ) - g(f) = e - V?) > 0. 

Therefore, dh(p,f) = g(f) and statement (c) follows. □ 



Fig. [2] illustrates Lemma 12.31 Because of Lemma 12.3} the travel time dh(p,f) between a 



a) 



b) 



Figure 2: a) If < a < ip v then all points move vertically to h. b) if ip v < a < j then some 
points move vertically to h, some other points move horizontally, and the rest of the points move 
directly to /. 
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demand point p and facility / simplifies to: 

' lb- /Ik, 



mm < 



IIP" 



Vlli + ^, 

Vlh i l|p "" /li2 
y 1 1 1 ' it 



(7) 



Given a solution (/, h) to the FFL-problem, we can always partition the set S of demand 
points into there sets S± := S±(f, h), S 2 ■= ^(/j h), and S3 := Ss(f, h) as follows. Set Si 
contains the points p G S such that x p < Xf and y p > yj and the points q £ S such that 
x q > xj and y g < y/. Set S2 contains the points p £ S such that x p < xj and p is either 
on or below /i, and the points q G S such that x q > xj and g is either on or above /i. Set 
S3 is equal to S \ (Si U S2). It is straightforward to obtain what follows (refer to Fig. [3]). 





a) 



Figure 3: Points p\, p 2 , and pz belong to Si, S2, and S3, respectively. In case a) we have 
< a < ip v . In case b) we have <p v < a < j. 



If < a < ip v then dh(p, f) is equal to: 



+ \^L^ll if p G Si 
if p G S 2 

Xf \ tana + „L^fJ ifpGS3 



1^ + |x P -a;/| tana , cosa , v 

\ x P~ x f\ 



\Vp-Vf\ - \x p -x f \t&na + c 



\y P -yf\ + \ x P- x f\^ ana + ^^ir "p^s 



and the objective function $(/, h) equals: 

w p\yp -yf\ + Yl w p\vp -yA-Y w p\ y p ~ y /i + 

peS-L pes 2 pes s 



X p - Xf \ + Wp|x p — x/| J + 



tana %>|£p — x/| — zu p 

— - — y~] w p\ x P - x f\ 

cos a • v L — ' 

Otherwise, if y>„ < a < j, then e^(p, /) is equal to: 

\x p - Xf\ + \y p - yj\ if pG Si 

|y P -2//l " l^-^/l^a + teF if PGS 2 
I K ~ x f \ - \y p - y f \ cot a + if p G S3 

and $(f,h) equals: 

^ tf p (\x p -x f \ + \y p -y f \)+^2 w p\Vp -Vf\ + ^2 w p\ x p ~ x ^ + 



P&S-L 



P65 2 

(r^^- tana ) E w p\ x p- x f\ + (—^- cota ) E w p\y P -yf\ ( 9 ) 



(8) 



7 



3 Reducing the search space 

Let G be the grid defined by the set of all axis-parallel lines passing through the elements 
of S. 




Figure 4: a) Lemma 1531 fa), b) Lemma [3711 (b) 



Lemma 3.1 There always exists an optimal solution (f,h) to the FFL-problem satisfying 
one of the next statements: (a) h contains a point of S and f is on a line of grid G, and 
(6) f is a vertex of grid G. Refer to Fig. [^} 

Proof. Let (/, h) be an optimal solution of the FFL-problem satisfying neither condition 
(a) nor condition (b). Using local linear perturbations, we will transform (/, h) into other 
optimal solution that satisfies at least one of these conditions. 

Assume first the case where ip v < a < ^ . Let fe > (resp. fe > 0) be the smallest value 
such that if we translate both / and h with vector (— <5i, 0) (resp. (62,0)) then / belongs 
to a vertical line of G or h contains a point of S. Given e E [— fe,fe], let fe an d h e be / 
and h translated with vector (e,0), respectively. Using Lemma 12.31 we partition S into 
four sets Z\, Z2, Z3, and Z$ as follows. Set Z\ (resp. Z2) contains the demand points 
doing a rightwards (resp. leftwards) movement to reach h. Set Z3 (resp. Z4) contains the 
demand points doing only a downwards (resp. upwards) movement to reach h. Observe 
that: 

e if p e Z\ 



dh E (p,fe) ~ dh(pJ) 



-e if p E Z 2 
e ■ c a if p G Z3 
— e ■ c a if p € Z4 



where c a = tana - co ^ v . Let Wi denote Y, P ez t w P' = 1, 2, 3, 4). Thus, for any 
e £ [— 61,62], the variation of objective function when we translate both / and h with 
vector (e, 0) is the following: 

$(f £ ,h £ )-$(f,h) = ^2w p -(d he ( P ,f £ )-d h (p,f)) 
= w p £ + w p(- £ ) + 

w P {ec a ) + ^ w p (-ec a ) 
p&z 3 P eZi 
= e(Wi- W 2 + c a (W 3 - W A )) 

Since (/, h) is optimal we must have 3>(/ e , h e ) — $(/, h) > for all e E [— <5i, fe]. It implies 
W1-W2 + c a {W 3 - W4) = and $(/ e , /i £ ) = *(/, fe) for all e E [-fe, fe]. Therefore, by 
translating both / and /i with vector either (— fe,0) or (fe,0), we ensure that / is on a 
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vertical line of G or h passes through a point of S, or both conditions. If it holds only that 
/ is on a vertical line of G, then we repeat the same operation in the vertical direction in 
order to ensure that / is on a horizontal line of G or h passes through a point of S, and 
condition (a) or condition (b) holds. Otherwise, if it holds only that h passes through a 
point of S, then it is straightforward to prove, by using similar arguments, that / can be 
translated along h in order to ensure that / belongs to a line of G, that is either vertical or 
horizontal, and condition (a) holds. In fact, for some demand point p of S, f will coincide 
with the point in which p enters h, that is, p' or p" . 

In the case where < a < <p v every demand point p moves vertically to h (Lemma I2.3|) . 
and we can proceed as follows by using arguments similar to the above ones. We first 
translate vertically both / and h with the same vector in order to ensure that h contains 
a point of S. After that, / is translated along h if necessary in order to / belongs to a 
vertical line of G and condition (a) holds. The lemma thus follows. □ 

Corollary 3.2 If a < ip v then there is an optimal solution satisfying Lemma HOI (a) . 

4 The algorithm to solve the FFL-problem 

Theorem 4.1 The FFL-problem can be solved in 0(n 3 ) time. 

Proof. We find an optimal solution by solving two cases separately. The first case is 
when solution satisfies Lemma 13.11 (a), and the second case is when solution satisfies 
Lemma (b) . 

In order to solve the first case we find for each demand point p and each line I of G an 
optimal angle a such that &(f a ,ha) is minimized, where h a is the line passing through 
p whose angle with respect to the positive direction of the x-axis is equal to a, and f a 
is the intersection point between I and h a . Assume w.l.o.g. that I is vertical and p is 
located to the left of I. It is easy to observe from equations ([8]) and Q that for any 
a E [0, j] the expression of $(/ a , h a ) has the form b\ + 62 tana + 63 cot a + — + 
where b±, b<i, 63, 64, 65 are constants. Furthermore, if we progressively increase the value of 
a from to j, that expression changes whenever sets Si, S2, and S3 change, that is, when 
h a crosses a demand point, f a crosses a horizontal line of G, or a = tp v . Then consider 
the sequence = ao < ai < ■ ■ ■ < a m = j of m + 1 = 0(n) angles, where each angle ojj 
(1 < i < m) is such that either h ai contains a demand point, f a . belongs to a horizontal 
line of grid G, or ctj = <p v . Notice then that the expression of $ (f a , h a ) is the same for all 
a G [ai, cui+i] (0 < i < m). If we preprocess the demand points S by constructing the dual 
arrangement of S [20], such a sequence can be obtained in 0(n) time by using both the 
Zone Theorem [20] and the order of the demand points with respect to the y-coordinate. 
Observe that if for a value of i we know the expression of &(f a , h a ) in the interval [ai, a^+i], 
then $(f a ,h a ) can be minimized in constant time in that interval. Furthermore, if h ai+1 
contains a demand point or f ai+1 belongs to a horizontal line of G, then the expression of 
^(/ai ha) in the interval [oj+i , aj+2] can be obtained in constant time from the expression 
of h a ) in the interval [a,, oji+i]. It is easy to see now that <J?(/ a , h a ), < a < \, can 
be minimized in 0(n) time by minimizing h a ) in [aj,Oi+i] for % = 0, 1, . . . , m — 1. 

Since there are n demand points and G has 2n lines, then an overall 0(n 3 )-time algorithm 
is obtained. 

We can proceed similarly in order to solve the second case. We find an optimal solution 
(u, h) for each vertex u of the grid G as follows. Let u be a vertex of G. Given an angle 
a, let h a be the line passing through u, whose angle with respect to the positive direction 
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of the x-axis is equal to a. Then, by Corollary 13.21 we look for an angle a E (<p v , j] 
such that the objective function <&(u,h a ) is minimized. It follows from equation Q that 
the expression of <&(u, h a ) has the form c\ + C2tana + C3 cot a + + for any 
a G (ip v , |], where ci, . . . ,05 are constants. If we progressively increase a from 93^ to j 
the expression of 3>(u, h a ) keeps unchanged as long as h a does not cross a demand point. 
The sorted sequence of values of a in which it happens can be obtained in linear time 
by using duality [20]. That sequence of values induces a partition of the interval (<p v , j] 
into intervals where in each of them the expression of 3>(u, h a ) is constant. We can now 
continue as was done above to solve the first case. Since G has 0(n 2 ) vertices then an 
overall 0(n 3 )-time algorithm is thus obtained. The result thus follows. □ 



5 A refinement of the algorithm 

Theorem 14.11 shows an algorithm that solves the FFL-problem by dividing the search of 
optimal solutions into two steps. It first looks in 0(n 3 ) time for an optimal solution that 
satisfies Lemma 13.11 (a), and after that looks within the same time complexity for an 
optimal solution satisfying Lemma l3.ll (b). In the following we show that for "reasonable" 
values of speed v we can simplify the algorithm of Theorem 14.11 by finding only optimal 
solutions that hold Lemma 13.11 (a) . We will use the next technical lemma. 

Lemma 5.1 Let a, b, c, and v > be non-negative constants and F : (0, — > M be a 
function so that 

1 — v sin x \ , / 1 — v cos x N 



F(x) = a + 6 : + c 

\ COSX / \ SIM / 

for all x £ (0, t;). Then next statements are true: 

(a) If a = and 6 = then F is constant. 

(b) If a > and 6 = then F is monotone decreasing. 

(c) If a = and 6 > then F is monotone increasing. 

(d) If a > and 6 > then F has no minima. 

Proof. Statement (a) is immediate. Let F' be the first derivative of F and observe that: 

„, . . / sin x — v \ , ( v — cos x 
F'(x) = a 5 +6 



cos 2 x J \ sin 2 x 

If a > and 6 = then lim x _ Sl o+ F(x) = a, lim £ - F(x) = —00, and equation F'(x) = 
has no solution in (0, ^) because v > 1. Therefore, F(x) is a monotone decreasing function 
and statement (b) thus holds. 

If a = and 6 > then lmx r _ 5 . + F{x) = —00, lim^x- F(x) = 6, and equation F'(x) = 
has no solution in (0, 5) because v > 1. Therefore, F{x) is a monotone increasing function 
and statement (c) thus holds. 

Consider a > and 6 > 0. Since linx r _ 5 .Q+ F(x) = lim^^Tr- F(x) = —00 it suffices to prove 
that equation F'(x) = has only one solution which must be a global maximum of F. 
Equation F'(x) = is equivalent to equation 

sin 2 x(v — sinx) 6 

G(x) := — r = - 

cos z x(v — cos x) a 
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We will prove that equation G'{x) = has no solution in (0, which implies that 
equation G(x) = | has a unique solution in (0, ^) because lim, !; „ !> o+ G{x) = < \ < 
+oo = lim„. 77- G(x). This will complete the proof. 

x ~ ' 2 

It is straightforward to see that equation G'(x) = is equivalent to equation H(x) := 0, 
where H(x) is equal to: 

2{y — sinx)(v — cosx) + sinxcosx(l — v(sinx + cosx)) 

Consider the function I : (0, 5) — > R so that 

I(x) = 2(v — sin x) (v — cos x) + sin x cos x(l — \/2v) 

for all x G (0, §). Since sinx + cosx < \/2 and sin x cosx > for all x G (0, 5), we have 
iJ(x) > /(x) for all x G (0, 5). We now show that the minimum of I(x) in (0, ^) is greater 
than zero, implying equation H{x) = has no solution. 

I'(x) = 2(— cos x(v — cosx) + sin x(v — sin x)) + 
(1 — v2v )(cos 2 x — sin 2 x) 
= 2(cosx — sin x) (cosx + sin x — v) + 
(1 — v2v) (cos x — sin x) (cos x + sin x) 

= (cos x — sin x) ( (3 — V2v) (sin x + cos x) — 2v 



If 3 — \[2v < then there is no x G (0, |) such that (3 — v / 2'w)(sinx + cosx) — 2v = 
because sinx + cosx is positive for all x G (0, 5). Suppose 3 — \/2w > 0. Then we have: 



(3 - \/2w) (sin x + cosx) - 2v < (3-V2v)V2-2v 

= 3V2 - Av 
< 

and thus there is no x G (0, |) such that (3 — \/2u)(sinx + cosx) — 2v = 0. Therefore, 
/'(x) = if and only if cosx — sinx = 0, that is, x = j. Let us prove that > 0. 




The roots of polynomial P(x) := 2x 2 — ^rpx + | are respectively equal to \ [^2^ — ^ 
and \ (^0- + X^) = Then we conclude that P(«) = I (f ) > because the main 



coefficient of P(x) is positive, is the greatest root of P(x), and v > ^j^. Since 
lmx E _ 5 . + /(x) = lim^Tr- I(x) = 2v(v — 1) > and I(x) > at the unique extreme point 
x = j, then I{x) > for all x G (0, f). Therefore, H(x) > for all x G (0, f), equation 
G'(x) has no solution in (0, ^), and then -F(x) has only one extreme point in (0, ^) which 
is a global maximum. The lemma follows. □ 

Lemma 5.2 If speed v is greater than ~ 1.060660172, then there always exists an 
optimal solution (f,h) to the FFL-problem satisfying Lemma \3.1\ (a), that is, h contains 
a point of S and f is on a line of grid G. 
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Proof. It suffices to show that for every solution (f,h), there exists another solution 
(/', h') satisfying Lemma IXTl (a) and <&(/', h') < <&(/, h). The proof is as follows. 

Let (/, h) be a solution of the FFL-problem so that h contains no demand point. Assume 
here angle a satisfies < a < ^. If a is such that < a < cp v , then the result follows 
from Corollary 13.21 Therefore, assume ip v < a < ^. We proceed to prove that there exists 
another solution (f',h') such that h! contains a demand point and <&(/', /i') < &(f,h). 
Consider the sets 52 = $2(/, h) and 53 = Ss(f, h) as defined before. 

Given an angle /3 € [(p v , \ ), let hp be the line passing through / so that the angle between 
hp and the positive direction of the x-axis is equal to [3. Observe then by equation ([9|) 
that $(/, h) is equal to: 

, , „ , , / 1 — v sin a \ , / 1 — v cos a \ 

$ /, K)=a[ +6 : + c 

\ cos a / \ sina J 

where a, b, and c are non-negative constants that depend on the coordinates of the demand 
points and the speed v. Then we can argue what follows by both noting that a > (resp. 
b > 0) if and only if S2 (resp. S3) is not empty and using Lemma 15.11 

We can rotate h with center / by either increasing or decreasing a to the value a' G 
[cp v , ^) in such a way solution (/, h a >) is obtained, where either a 1 = (p v or h a > contains a 
demand point of S2 U S3, and < ^(/, /i). If a' = fv then the result follows from 

Corollary 13.21 Otherwise, if h a i contains a demand point of S^US^, then (/', h') = (/, h a /) 
is the desired solution, and to finalize the proof, we translate /' along h' , if necessary, as 
was done in the proof of Lemma l3.1| in order to ensure /' belongs to a line of grid G. The 
result thus follows. □ 



6 Examples 

In Fig. Oand Fig. Owe show the same example, consisting of nine demand points p\ , P2 , • • • , P9 
Each demand point is represented by a solid dot, and labeled with a triple, the first two 
components are the coordinates and the third component is its weight. Facility point 
/ is represented by a cross. In both examples optimal solutions satisfy Lemma 13.11 (a), 
that is, highway contains a demand point and facility point belongs to a line of grid G. 
As expected, when we increase the highway's speed from the example in Fig. [5] to the 
one in Fig. [6l the shortest paths to the facility point change according to the claims of 
Lemma 12.31 



p s = (S,5,8) 




Figure 5: The highway contains a demand point and facility point belongs to a line of grid G. 
Some points perform an horizontal movement to reach the highway. 
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58=^(8,5,8) 




Figure 6: The highway contains a demand point and facility point belongs to a line of grid G. All 
points perform only a vertical movement to reach the highway. 

In Fig. [5] speed v is equal to 1.2 and highway h of the optimal solution contains point 
Ps, facility point / = (— 1, —0.14335) is on the vertical line passing through p§, and there 
are some demand points moving horizontally to reach facility point /. The value of the 
objective function <£(/, h) is equal to 525.83. 

In Fig. [6] speed v is equal to 1.5, highway h of the optimal solution contains point p2, 
facility / = (—1,-0.7758) is on the vertical line containing p§, and all demand points 
move vertically only to reach the highway. Since speed is greater than speed in Fig. El the 
value of the objective function 3>(/, h) reduces to 471.55. 

In Fig. [7] we present a different example consisting of nine demand points p\,...,pg, with 
the aim of showing the existence of configurations for which optimal solutions satisfy only 
Lemma [3. II (b). In this example speed v is equal to 1.04 < ^j^. Highway h of the optimal 
solution contains no demand point and facility / = (0, 0) is on a vertex of grid G, in fact, 
it is located on both the horizontal line through pi and the vertical line through pg. The 
value of <£(/, h) is equal to 336.2. 



pi = (- 




Figure 7: The highway contains no demand point and facility point is a vertex of grid G. 

7 Conclusions and further research 

We have solved in 0(n 3 ) time the problem of locating at the same time a facility point and 
a freeway of variable length, among a set of demand points, in order to minimize the total 
weighted travel time from the demand points to the facility. Some examples are presented 
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to show that there exist optimal solutions corresponding to each type of solutions that our 
algorithm considers. 

A natural restriction to be considered in further research of this problem is to upper 
bound the length of the highway, that is, to consider that highway has fixed length. In 
this case, there also exist optimal solutions in which the facility point belongs to the 
highway. This was in fact showed in Proposition 12.21 because in the proof we did not 
change the lenght of the highway. It is not hard to see that when the highway's length is 
fixed, the shortest paths from the demand points to the facility point can be discretized 
as follows. If < a < ip v , then we distinguish three regions R\, R2, and i?3 as depicted 
in Fig. [8] a). Points belonging to R\ U i?3 move to the nearest endpoint of h, and points 
of R2 move vertically to h. Otherwise, if ip v < a < ^, then eight regions R\,. . . ,R$ can 
be identified as shown in Fig. [8]b). Points of R± U R2 move to the nearest endpoint of h, 
points of R3 U R4 move directly to /, points of R§ U Rq move horizontally to h, and points 
of Rj U i?8 move vertically to h. 



a) 




Figure 8: Discretization when the highway's length is fixed. 



We believe that with the above discretization, the search space of optimal solutions can 
be simplified by using a similar (and more detailed) statement as Lemma 13.11 This will 
permit to obtain an algorithm similar to the one presented in Theorem 14.11 

Other variant to be considered in further research is the problem of locating at the same 
time the facility point and a turnpike of variable length. This will extend [12^14] , 
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